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Abstract 



^ . Motivated by the renormalization group (RG) approach to c = matrix model of 



Brezin and Zinn-Justin, we develop a RG scheme for c = 1 matrix model on a circle 
and analyze how the two coupling constants in double scaling limit with critical exponent 
flow with the change in length scale. The RG flow equations produce a non-trivial fixed 
point with the correct string susceptibility exponent and the expected logarithmic scaling 
violation of the c = 1 theory. The change of world-sheet free energy with length scale 
indicates a sign change as we increase the temperature, indicating a phase transition due 
to liberation of the non-singlet states. At low temperature, the RG analysis also lead to 
T-duality of the singlet sector free energy. The RG flow to the c = 1 fixed point can be 
understood as the decay of unstable DO-branes with open string rolling tachyon to the 2D 
closed string theory described by the end point of the flow. The amplitude of the decay 
is extracted from the change of the world-sheet free energy described by the RG process 
and is in accordance with the prediction from boundary Liouville theory. 
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1 Introduction 

Back in the late 80's and early 90's, there had been extensive work on low dimensional string 
theories, namely the c < 1 non-critical (or D < 2 critical) bosonic string theories, both in 
the discretized large N matrix model approach |], |2|, |3], |J and in continuum Liouville theory 
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approach 0, 0] by looking at them as 2D quantum gravity coupled to c < 1 matter. One 
of the motivations was to learn something about the non-perturbative effects of 2D quantum 
gravity and string theory and to apply that knowledge to higher dimensions. Studying 2D 
gravity with simple matter was thought to be useful to probe issues of topology change in 
a finite and tame gravitational theory and to construct a string theoretical basis of QCD. A 
remarkable fact is that all the ideas of developing non-perturbative techniques using branes and 
dualities, string field theory etc. is now going back to matrix model. According to the recent 
idea, initiated in ||, the eigenvalues of c = 1 matrix model represent unstable D-branes in the 
dual two-dimensional string theory. Before going into the details of this idea in the context of 
this paper, we describe briefly some of the difficulties in understanding c = 1 matrix model on 
compact space to sketch the main motivation of this paper. 

In the discrete approach (see |], |IU], |TT], [12], I3[ [13]] for review), the Euclidean path-integral 



for two dimensional gravity can be written as a sum over discretized random surfaces and is 
evaluated using the large N or the planar limit of the appropriate random matrix models. 
However the large N limit, the lowest order in the string perturbation theory, is not the only 
thing that we can extract from the expansion of the world-sheet free energy. In the continuum 
limit, for any genus, the area of the surface is large compared to the elementary polygons of the 
discretized surface. Hence for any fixed genus G, i.e., at order N~ 2G in the 1/A-expansion, we 
should stay in the vicinity of the singular point g = g c of the free energy at which the mean area 
diverges. Thus the limit of interest involves large N and small (g — g c ) with g 2 s = Ar2( . g _ g 1 ^ 2 _ 7str 



(where 7 sir is the string susceptibility exponent) fixed. This double scaling limit | jr5| , |TBj , |TT 
led to the solution of the matrix models to all orders in string perturbation theory. The genus 
expansion of the world-sheet free energy can be written as 

oo 

F(N,g) = J2N 2 - 2G (g-g c Y 1 ~ G ^^. (1.1) 

G=0 

Also the genus expansion becomes universal, i.e. almost entirely independent of how the surface 
is discretized. 

The simplest model involving integral over one matrix variable has been solved to describe 



pure two-dimensional gravity (the c = case) [|15|, [16], |17|]- The one-dimensional hermitian 
matrix chain models have been solved in the double scaling limit and identified with the c < 1 
minimal models coupled to two-dimensional gravity [|Tj| [HJ |2(], 2l|. For this general case of (p, q) 



minimal models coupled gravity, Douglas has proposed a solution in terms of the generalized 
KdV equations , where the non-perturbative partition function is given by the square of the 
r-function of the KdV hierarchy, satisfying the string equation. The hermitian matrix quantum 



mechanics, describing c = 1 theory, has been solved p2, 23| £24, and is interpreted in terms 



of two-dimensional string theory p6l p7] , where the role of the extra dimension is played by the 
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conformal factor of the world-sheet quantum gravity. 

The solvability of the c < 1 matrix models of N 2 degrees of freedom is mainly due to the 
fact that only N eigenvalues contribute, and the other N(N — 1) angular degrees of freedom 
decouple. The SU (N) singlet spectrum is described by N free fermions moving in a potential. 
The non-singlet states, coming from the angular degrees of freedom have energy diverging as 
logarithm of the cut-off. So they are pushed to infinity in the continuum limit. As a result, the 
singlet wave functions of the c = 1 matrix model are sufficient to account for all physical states 
of the continuum string theory in a non-compact dimension, or in a circle bigger than a critical 
size, R > R c ||28|| . This decoupling corresponds to the confinement of Berezinskii-Kosterlitz- 
Thouless (BKT) vortices |29], |3(J , which wind around the target space plaquette of lattice size. 
Although each of the vortices is suppressed by a divergent action in the continuum limit, its 
entropy factor, given by the number of places it can be found, grows and for sufficiently small 
circle R < R c , a Kosterlitz-Thouless phase transition occurs. As a result, although exact results 
in matrix model have been proved to be very powerful in describing low dimensional string 
theory in uncompactified dimension, understanding physics from very small compact target 
space is lacking. Although, for c = 1 case, X and the extra hidden dimension or the Liouville 
field are a priori spatial dimensions, it is more convenient to continue X — > iX, and denoting 
t = iX as the Euclidean time. If we have t wrapped on a circle, t ~ t + 2iTnR, we have c = 1 
string at finite temperature. 

For c > 1, all the iV 2 angular degrees of freedom become relevant and the situation becomes 
difficult to pursue in the matrix model framework, until we understand considerably how to 
deal with the non-singlet states. At this moment the predictions of continuum theories for the 
well- understood c < 1 cases turn out to be meaningless for c > 1, as the KPZ-DDK formula 
H 01 f° r the string susceptibility exponent and for scaling dimensions of matter operators 
lead to nonphysical complex values. Having said all these, the difficulties with c > 1 models are 
not related to just these technical issues. The continuum approach indicates that these models 
are tachyonic [[H], Light-cone quantization of certain c = 2 matrix model gave some partial 
results p3| , but at the same time faces potential difficulties as the non-singlet states do not 
confine even if non-compact target space is considered. 

The motivation of this paper is to initiate an approach to deal with the less understood non- 
singlet sector in matrix quantum mechanics on a compact space. We develop a renormalization 
group (RG) scheme, generalizing the work of Brezin and Zinn- Justin |34]] on c = to the c = 1 
model on a circle. The motivation is to develop a scheme which would reasonably reproduce the 
known cases that are solvable and to understand, at least qualitatively, the physically interesting 
situations which can not be simplified because of the lack of the solvable structure. The RG 
analysis of the c = model, the simplest of the solvable examples, studied in has been 
extended to various general cases (see for example []35], |36|, [37|, ^]). In |36|| , reparametrization 
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invariance and the the loop equations are used to eliminate some of the induced interactions 
which appear in the RG transformation at higher orders in coupling. However extension of this 
method to matrix quantum mechanics on circle is not obvious. We instead directly generalize 
[|34|j by explicitly evaluating the determinant obtained by integrating out part of the matrices 
to get a RG flow in the space of actions in general representation. This enables us to study the 
models which do not allow angular integration, such as the c = 1 model on small circle, or c > 1 
matter coupled to gravity. 

As an illuminating interpretation and feedback from our RG analysis, we will digress to the 
recent progress in understanding the dynamics of the boundaries in the 2D string theory, that 
realizes the duality between c = 1 model and 2D quantum gravity coupled to c = 1 matter 



as an exact open/closed string duality ||, |39|, [40], fy], [42], [43], [44]]. According to this idea, the 
quantum mechanics of SU(N) invariant matrix variables in an inverted oscillator potential is 
visualized as the quantum mechanics of open string tachyons attached to N unstable DO branes 
that decays into (i.e. dual to) Liouville theory coupled to c = 1 matter describing 2D closed 
string theory together with its DO branes. In matrix model, an unstable brane corresponds to 
a free fermion excited to the top of the potential (in presence of the rest forming the static 
Fermi sea). It then decays to the closed string vacuum by rolling down to the Fermi level as an 
unstable eigenvalue trajectory 

z{t) = y/2iwt! A e* . 

The amplitude of this decay calculated by continuum methods can be read from the matrix 
model analysis considering the bosonization of the relativistic fermions in the asymptotic limit 
pCfl. Alternatively as proposed in JO], one can formally view the Vandermonde determinant 



corresponding to the (N + l)-th unstable eigenvalues in presence of the rest N, to be related to 
the bosonization field and thus to lead to the decay amplitude. 

To understand from our RG analysis, we consider integrating out one of the eigenvalues in 
presence of the static Fermi sea of the rest of the N eigenvalues. The decay amplitude is shown 
to be contained in the change of the world-sheet free energy given by the ratio Zn + i/Z^ of 
the partition functions. The nice thing is that the end point of such a decay is explicitly a 
c = 1 fixed point of the flow. We compute the above-mentioned Vandermonde determinant, 
the entropy arising out of the loss of information due to integrating out one of the eigenvalues, 
taking the unstable eigenvalue trajectory to be that corresponding to the ZZ-boundary state 
4^j tensored with Sen's boundary state JKJ. This leads to the explicit relation between the two 



different views in [|I§ and (41| to compute decay amplitude and reproduces the expected time 
of decay 

At ~ In A . 

Understanding the dynamics of the boundaries in the 2D string theory, i.e. the D-objects 
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of the boundary Liouville theory, from the RG analysis of general (N x N) matrices on a 
circle is much more interesting. In this case, by tuning the matrix coupling constant and the 
mass parameter with N — > oo in the double scaling limit, one tunes the bulk and boundary 
cosmo logical constants respectively to arrive at a relation between the two. This is given by 
the integration of the flow hitting the c = 1 fixed point. Such a relation presumably represent 
Dirichlet or Neumann boundaries that are present in the Liouville theory coupled to c = 1 
matter. We will show that for large circle we arrive at such trajectories with correct scaling 
between the two cosmological constants. However a detail analysis of such trajectories for large 
and small circle will be extremely interesting and will be studied in detail in a future publication. 

The organization of the paper is as follows. In section 2, we present the details of the 
RG calculation of matrix quantum mechanics on a circle with a cubic potential and obtain 
the set of beta function equations for the couplings. In section 3, we analyze the flow, the 
fixed points and the critical exponents. As understood previously, the usual KPZ-DDK scaling 
laws [§, H [7| for c < 1 are obeyed with a logarithmic scaling violation for the c = 1 model 
23, 1231, EM, 1231, ESI. With the understanding of the RG flows, here we are able to obatin the 



correct critical exponents of the c = 1 strings, and to reproduce the scaling behavior for the 
free energy with expected scaling violation for the singlet sector. Also the RG trajectory near 
the nontrivial fixed point shows the correct scaling behavior between the bulk and boundary 
cosmological constants. From the running of the prefactor of the partition function written in 
the renormalized couplings, analogous to the running due to the wave function renormalization, 
the free energy is observed to change sign near R = 1 for small value of the critical coupling. This 
is reminiscent of the Kosterlitz-Thouless transition at self-dual radius triggered by the liberation 
of the world-sheet vortices. In section 4, we discuss that in what sense the RG analysis can 
capture the T-duality respected by the singlet partition function. In section 5, we discuss how 
the change of the world-sheet free energy in integrating out one matrix eigenvalue can be seen 
in the recent context to lead to the amplitude of the closed string emission from the decaying 
brane. In section 6, we conclude with some open questions. 



2 The Basic Set Up and the World-sheet RG Calculation 

The existence of double scaling limit indicates that a change in length scale induces flow in the 
coupling constants of the theory in a way that one reaches the continuum limit with desired 
critical exponents. In this continuum limit, as the matrix coupling constant approaches a critical 
value g c the average number of triangles (no) in triangulations at any genus G diverges with the 
exponent —1 while the length of the triangles or the regularized spacing of the random lattice a 
is scaled to zero by taking iV — > oo to keep the physical area a 2 (no) or equivalently the string 
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coupling g s fixed, 



9^ 9c (n G )^(l-G)( l0 -2)(l-g/g c ) 



-i 



oo , 



N -> oo =>- a ~ N 2 -^o -> , 
with a 2 (n G ) ~ JV -3 =ro (1 - g/g c ) = const, or, g~ 2 = N 2 (g - g c ) 2 ~ 10 = const . (2.1) 

Thus one would naturally try to understand how the two parameters of the theory, the size 
of the matrices N and the cosmological constant (mapped into the matrix coupling g), evolve 
at the constant long distance physics with the rescaling of the regularization length in the 
triangulation of the world-sheet. The flow equations will automatically give rise to the correct 
scaling laws and the critical exponents around the nontrivial IR fixed points governing the 
continuum physics. In the Wilsonian sense this can be achieved by changing N —> N + 5N 
by integrating out some of the matrix elements, which is like integrating over the momentum 
shell A — dA < \p\ < A, and compensating it by enlarging the space of the coupling constants 
g — > g + 5g. Here the space of coupling constants will contain both the matrix coupling g and 
the mass parameter M 2 . 

Following the RG scheme of Brezin and Zinn- Justin we construct the flow equations by 
integrating out a column and a row of an (N + 1) x (N + 1) matrix, reducing it to an x N 
matrix. One expects the process to lead to the following key relation satisfied by the matrix 
partition function 



Z N+1 (g, M, R) = [\(g, M, R)fz N tf, M', R!) 



(2.2) 



with 



M 



9 



12 




(2.3) 



and 




(2.4) 



Then the string partition function 



g, M, R) — -^r In Z N (g, M, R) 



(2.5) 



satisfies the Callan-Symanzik equation 
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f) 

N— - (3(g, M, R)— - (3(g, M, R)— + 7 (p, M, R) FiX.g. M. R) ■■ rfjj.M.B). (2.0) 



2.1 Integrating out a row and a column 

Let us consider the (N + 1) x (N + 1) matrices <pN+i(t) and decompose them into N x N 
matrices <f>N(t), iV-vectors v a (t) and v*(t) (a = 1, . . . , N) and a scalar a. For the time being we 
can choose a = as they are of relative order 1/N and can be ignored in the double scaling 
limit. 



<!>N+i(t) 



4>N(t) V a (t) 



(2.7) 



V*a(t) 

For simplicity, we consider a cubic potential, which will serve as a wall stabilizing the inverted 
oscillator potential. In terms of iV + 1 dimensional matrix variable, the action reads 



S N+1 [<j> N+1 {t),g,M,R] = (iV+1) 



2ttR 



dt Tr 



l^N + i(t) + lM 2 ^ 2 N+1 (t)- 9 -^ N+1 (t) . (2.8) 



The parametrization ( |2.7|) gives the simple relations, 



Tr^ +1 = Tr0^ + 2k v*^~ z v + 0(v 4 ) , 
Tr^ 1 = TY<p k N + (2k + 1) v*<f^- x v + 0(v^ 2k - 3 " 



(2.9) 



The higher order terms in v*v can be neglected as they are supressed by powers of 0(1/ N). 
The resulting partition function can be written as 



Z N+1 [g,M, R] 



V N2 (f) N (t) e-W+ifbf?'** {|^(*)+|M 2 ^(*)-f^(*)} 

<f> N (2TrR)=4> N (0) 

x f V N v(t)V N V*(t) e -( N + 1 ) Io WH dt v*(t){-d?+Ml-g<t> N (t)]v(t) 

Jv,v*(2irR)=v,v* (0) 



(2.10) 

The above partition function is identical to the one considered in |fl?| , where the c = 1 
matrix model suitable for Veneziano type QCD has been considered to study open strings. 
Both color iV and fermion (quark) flavor Nf has been taken to be large. Here the quarks are 
bosonic (the vectors). It is precisely these fields in the fundamental representation of the global 



S 



SU (N) group, which generate boundary terms in the Feynman diagrams. As usual, we will use 
adiabatic treatment of first integrating out the quark loops. Integrating over the quarks, we get 



Z N+1 [g,M,R] = (^—) N [ V N2 Mt) 

WV+l/ J^ N (2^=^(0) 

xexp{-JJ dt [(N + l)Trji^(t) + l -M 2 <f N (t) - 9 -<f N (t) 



+N Tr log{-^ 2 + M 2 — g<f> N {t)}] } . 



(2.11) 



Logarithm with minus sign arise if the integration is performed on iV flavors of fermions. 
Although such model has been considered as matrix model for open strings, where the logarith- 
mic term has the effect of generating boundaries in the world-sheet, there are some differences. 



Unlike the zero dimensional case [^, |50], |51|| , this logarithmic term does not arise by integrating 
out fields (quarks), which are lxN matrices that couple to 4>n- As a result, in those open string 
models, the couplings in front of the logarithm and in its argument are introduced by hand, 
rather than being determined by the original closed string action. One tunes the couplings to 
some appropriate values and the adiabatic treatment of first integrating out the quark loops rise 
to tearing phenomena. Interesting critical behavior of getting phases with torn surface in the 
c = case j|9|| also occur in the one- dimensional case [|7j if the dynamical loops are generated 



by bosons without kinetic term. Absence of the kinetic term makes fermion loops uncorrelated 
in time. One can ignore the derivative term inside the logarithm if the mass and the couplings 
are large enough. But there are interesting critical behavior when the argument of the logarithm 
without the kinetic term approaches zero [^, [17]. In jPJ, c = 1 model with explicit expression 



for the fundamental fields has been considered and is the one which is closest to our model 
obtained after integrating out one row and one column of flavor degrees of freedom. If one 
considers infinite line, only ground state is relevant. For some choices of the coupling and the 
mass of the particle moving around the boundary of the holes, it has been possible to find the 
ground state |48|| , but the exact spectrum is not known. We will return to this discussion in 



section 5, when we will interpret the results from the RG analysis in the recent context of time 
dependent tachyonic decay of D-branes in two-dimensional string theory. 

Now we return to evaluate the determinant by standard Feynman expansion. In general 
it is a non-local object. For simplicity one could consider the constant <pN mode, which is 
equivalent to studying the effective potential to determine the phase structure. But here we will 
treat more general case of evaluating the determinant with flavors coupling to general <$n by 
performing the calculation in Fourier transformed variables with discrete momenta arising from 
the compact target space. The complicated induced interactions arising from the logarithmic 
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term are ignored in small field approximation in order to get back the beta functions that depend 
only on the original couplings. 

Rescaling the vectors v(t) — > — ^ v ^ N+ ^ > the 11 dependent part of the partition function turns 
out to be 



I[g,M,<f>N,R] = WWTT^ [ V N v(t)V N v*(t) 

r r 2nR dt i 
exp [-J — + ^*(M 2 - ^)i;}J , (2.12) 

Before integrating out the vectors to get the determinant det 1 1 —d 2 +M 2 —gcp^l | _1 and expanding 
it in Feynman diagrams, it is convenient to Fourier transform all the fields as 

0(t) = V O m O m = -^e^t), (2.13) 

m ^oo ■><> ^R 

with 



5 mn = / e « , 5(t-t) = > e n yt 1 

J 2nR ' V ; 2nR ^ 

J[} m=-oo 

In terms of the Fourier modes, the f-integration can be expressed as 



I[g,M,(b N ,R] = 1 = — [ (T\dv*dv*) e -Zm" m m {w+ M2 ) v <"+9Z m ,i< 

1 J y2itR(N + l)Y N J VA n A n nl 

(2.14) 

where we have neglected the 0(1/N) terms. 
2.2 One loop Feynman Diagrams 

In order to carry out the v integration diagrammatically, let us now define the following operators 



°tl = (^ + M 2 )5 mn , 0^1 l {gA)=g<t>m-i, (2.15) 

The inverse of these operators define various propagators and vertices according to figure |l|. 
Hence the integral becomes 

I[g, M, fa, R, N] = —-L— exp [ - £ 0%_ h (g, 0)] I (R 7 M, N) , (2.16) 
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(b) 



Figure 1: The propagators and vertices: (a) [C^n] S O 3 ) l^m-iidi 



where the gaussian part is as follows 



I [R,M,N] = J (Jldv^dvj) exp[- Y,<°%&»] 



[2.17) 



In order to perform the gaussian integration, we rescale v n as v n — > v n /(n 2 /R? + M 2 ) 1 / 2 and 



use (|2.13|) . Performing the gaussian integration, we get 

1 ' J Vsinh 2 7ri?M/ 11 \nJ 
where we have used the standard relation 



(2.18) 



n=l 



7/2 



TTX 



no v 

- LJ - V n z / smh7rx 

n=l 

Inserting this into ( [2.16|) , the v- integration becomes 



I[g, M, <t> N , R, N] = C(R, M, N) Efo, M, 0^, i2, iV] 



where 



C(i2, M, iV) 



7T 



3 R 2 



N oo 



n=l 



n<-> 



(2.19) 



;2.20) 



.(iV + l)sinh 2 nRM_ 

In ( 2.19|) , M, 07v, i?, TV] represents sum of one loop Feynman diagrams as shown in figure [| 
The sum T,[g, M, <f> N , R, N] can be expressed as 
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Figure 2: The diagrams contributing to £ at one-loop order. 



E[g,<f> N ,R,N]=l+gTi 

1 



E^ 



L n R 2 



+ M 2 



+ ^Tr 
2 



E 



y m—nr"n—m 



g_ 

3! 



Tr 



E 



- (^ + AP)(^ + Af 2 )( ^ + M 2) 



<p m -l<pl- n (/> rt 



0(g 4 ) (2.2i; 



2.3 Evaluation of the diagrams 



In order to evaluate the one loop correction to the effective action, we inverse transform the 
Fourier modes according to the rule ( |2.13|) and sum-up the set of infinite series using the formulae 



sum r 



exp[i(m/R)t] itR cosh(vrM J R - Mt) 



M 



sinh it MR 



< t < 2-kR . 



[2.22) 



and, 



E 



m/R exp[i(m/ R)t] mR smh(iraR — at) 



sinh TiaR 



iirR , 
—itxR . 
0. 



< t < 2nR. 

t -> 0, 
t -> 2ttR, 
t = = 2nR . 



(2.23) 



These summations give the corrections to the coefficients of the various terms in the ac- 
tion, after E[g, M, (j) N , R, N] is exponentiated and log-expanded using small field approximation. 
Since after performing the inverse Fourier transform, the various terms has nonlocal integrals 
over several one dimensional dummy time variables, we breakup the variables into center of mass 
and relative coordinates. Then we expand the functions about the center of mass coordinates, 
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assuming the relative coordinates to be small enough, and consider integration over the relative 
coordinates. 

The details of the diagram evaluation are given in the appendix taking into account all 
the above considerations and evaluating the integrations over center of mass and relative time 
variables. Collecting all the contributions up to 0(000) from the appendix, the expression for 
£[(?, M, 0w, R, N] becomes 



p2wR p2wR 

%, M, <j> N , R,N] ~ 1 + F gl (R, M)g dt Tr</> N (t) + F g2 (R, M) g 2 / dt Tr0^(t) 

Jo Jo 



/2-nR i i-2-kR i 

dt -Tr0^(t) + F g3 (R, M) g 3 J dt -Tr<j> 3 N (t) . 



Note that, in the evaluation of S we keep the contribution from the nonlocal terms of the action 
up to the kinetic term and the contribution from the higher order terms in the matrix field up 
to the cubic term. All other operators are redundant for our purpose and are negligible due to 
the small field approximation. 

In the above expression, F(R)s are defined as follows 



F g i(R, M) = ^cothTiMR, 

F q2 {R,M) = -= (-7TMRcosh27rMR + -smhA7rMR) , 

1 (TlMR 7T 3 M 3 i? 3 

F q2 (R, M) = ^ ( cosh AttMR cosh lixMR 

9 K ' M 5 sinh irMR V 16 6 

- — (1 + 8n 2 M 2 R 2 ) sinh 4nMR^ , 

,. ttMR 
F 9 s(R,M) = ,^,,^3 



64M 5 sinh 3 7rMi?(cosh 2ttMR + cosh47rM J R) 
[An MR (3 cosh irMR + 2 cosh 3ttMR + 2 cosh hnMR + cosh 7rrMR) 

+ sinh irMR + sinh 3nMR + sinh 5nMR + 2 sinh 7ttMR + sinh 9irMR] . 

(2.24) 



2.4 Elimination of the tadpole term 

The term proportional to Jq 7tR dt 0jv(t) is unwanted. As usual, in order to remove this term, 
we change the background 0at(£) — > 0jv(t) + /, and set the net coefficient of the term linear in 
to zero. This fixes the value of / as 
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J 2V N N ) l\ N N J U N N J 

and accordingly modifies the coefficients of all the terms in the action. After accommodating 
all the changes, the expression for Z N+ i turns out to be 



J N+1 



j V N2 (p N exp 


p2nR 


- N Tr J dt j 







9 2 F 9 2 \ <j> 2 ( M 2 g 2 F g2 g 2 M 2 F gl x 2 (t) 
JV / 2 I JV JV AT 72 



ff 3 F g3 \0 3 (t)) 
iV / 3 J. 



1 + 



(2.26) 



where the expression for J^"(<?, M, R, N) is given by 



&(g, R,N) = 9 -^f+(M* + — + °-**) i -+(g + ^- + 9 -^n) f - , (2.27) 
vy '' ; N \ N NJ2 V N N J 3 V ; 

which is of 0(^2 ). 



2.5 Rescaling of the fields and the variables 

In order to restore the original cut-off, we perform the following rescalings so that the effective 
action is of the same form as the bare one but with the renormalized strength of the couplings. 

<j> N (t) -> p<p' N (t') , t' -> t(l - h dl) , R 1 -> R(l - h dl) , (2.28) 

where dl = 1/N, and set the overall coefficient of the kinetic term to zero. Here we can assume 
that 

h = ^*jg i M i h i j(R). (2.29) 
hi 

The functional form of h can be guessed from the contribution of the Feynman diagrams in the 
behavior of the flow near the fixed points. This sets the value of p to 

p=l + ^(h-l + g 2 F g2 )dl + 0(dl 2 ), . (2.30) 
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2.6 The beta function equations 

The effective action is of the same form as the bare one, but with a renormalized strength of 
coupling. The resulting partition function is given by 



Z N+l = X' N2 J V N U' N (t') exp[-ArTr^'dt'(^^ + M' 2 ^^-^^ 3(t/) 



3 



(2.31) 

where, 

A' N * = C{R, M, N) exp[-2vr RN 2 ^{g, M, N, R)]p N2 . (2.32) 

Neglecting the 0(dl 2 ) terms, the renormalized coupling and mass are expressed in terms of the 
bare quantities as follows: 



9 

M' 2 
A' 



g+(^h-^jgdl+ [F g3 (R, M) - ^F g2 (R, M)]g 3 dl , 

M 2 + [2hM 2 + g 2 (l- M 2 )F g2 (R, M) - g 2 M 2 F g2 (R, M)]dl . 
tc 3 R 2 l „ 1 



In 



sinh ttRM 



dl + -[(h-l)+g 2 F g2 (R,M)]dl. 



[2.33) 



Here, in simplifying the part C(R, M, N) n in the expression for A', we have assumed that for 
any value of R, 

r 3 D 2 _ r ^fftn 



C(R, M,N)- 



exp 



1 + 



N 11 V sinh it RM ) 



(N + l)(n\ 



n— >oo, N—kx 



nH^rm) +0 ^- 

Also, the term exp[-2irRN 2 ^(g, M, N, R)} contributes only a factor of 1 as JP{g, M, N, R) 
0(dl 2 ). The resulting beta function equations are given by 



7T 3 R 2 



(2.34) 



P 9 



dl ~ 

dM 2 

~dT 

d\ 

~dJ ~ 



2hM 2 + g 2 [(l - M 2 )F g2 (R, M) - M 2 F g2 (R, M)) 



In 



7T 



3 R 2 



sinh ttRMI 



[(h-l)+g 2 F g2 (R,M)]. 



The relation (|2.28 ) indicates in some sense a renormalization of the radius R as 



R dR hT7 

P R = - = -hR. 



[2.35) 



[2.36) 
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Now before going to the detail analysis of the fixed points let us try to understand few 
things about the structure of the beta function equations. The much of the structure depends on 
understanding the quantity h. One can clearly see for g = the gaussian model is never expected 
to flow and hence for such a fixed point h = (i.e. corresponding to the trivial rescaling t' = t 
and R' = R). The situation is different for a non- vanishing h. Demanding the mass parameter 
to be always at some fixed value of M (we will use M 2 = 1 for simplicity) and (3 m 2 = 0, one can 
easily determine some h = h(R) for nontrivial fixed points g* ^ 0, as there are two independent 
equations (f3 g = 0, (3m 2 — 0) and two unknowns (h(R) and g). This is extremely interesting as 
it could indicate a phase transition at certain radius due to turning on (i.e. being relevant) the 
operator coupled to the cosmological constant. Keeping M 2 = 1 for simplicity is consistent with 
the value of the mass parameter (M 2 = ^7) one originally works with in the matrix partition 
function to visualize the matrix path integral as the generator of the discretized version of the 
Polyakov path integral of 2D bosonic string. In recent identification of the matrix quantum 
mechanics with the quantum mechanics of open string tachyon on unstable DO-branes the mass 
parameter M 2 = ^7 is identified with the open string tachyon mass. However, we will use a 
framework of the flow of a general M 2 (i.e. h = )to discuss the presence of the boundaries. 

Another point to be mentioned is that in this general analysis though we will see the signature 
of the singlet sector to be more apparent, the non-singlet sector does render indirect signature 
on the flow. This is reflected in the R dependence of the change of the world sheet free energy 
with the scale, i.e. in the (3\. For very small g* it exhibits a sign change at certain radius 
which is reminiscent of the tendency of the non-singlet sector to be liberated at certain critical 
temperature (R c =1). To capture this effect one needs to turn on the right operator coupled 
to the fugacity of the vortices. In a separate analysis we will work with a gauged matrix model 
with a gauge breaking term that couples to the fugacity of vortices and drives the system to 
the phase where the vortices are liberated, above a critical temperature. 



3 The c = 1 Fixed point 

We will now analyze the fixed points of the flow equations given by the simultaneous solutions 
of P g = (3 M 2 = 0. As we discussed before, for the trivial rescaling of the coordinates and the 
momenta, t' = t and 1/R' = 1/R , one can assume h = 0. In this case, the gaussian fixed point 

= (0, 0) satisfies both the equations (3 g = and (3 M 2 = trivially. The nontrivial fixed point 

is given by, 

g* = ±J — ; , (3.1) 

V 2F g3 (R,M*)-3F g2 (R,M*) 
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where M* is determined by the equation, 

g* 2 ((l — M* 2 )F g2 (R, M*) — M* 2 F g2 (R, M*)) =0. (3.2) 

We will see in the next section that this nontrivial value of g* will always characterize c = 1 
fixed point as it gives the critical exponent of c = 1 for any R for h = 0. Using the value of g* 
at the nontrivial fixed point Q3.1| ) in the equation ( |3.2| ), one has 

(1 - M* 2 )F g2 (R, M*) - M* 2 F g2 (R, M*) =Q 
2F g3 (R,M*)-3F g2 (R,M*) 

For any value of R, this equation trivially has solution for very small value of M*. However, for 
large R there are also solutions for large M*. In the large R limit the equation (|3.3| ) takes the 
form 

^/8n 2 M* 2 R 2 -8M* 2 -4nM*R + 9\ , s 

M I = (3.4) 

V 24vr 2 M* 2 i? 2 - 10ttM*R + 3 J 1 ; 

The solution for large M* corresponds to the situation that the denominator is much larger 
than the numerator (which solves with a precision of 70% ). This gives 

M* > m+ , or M* < m_ , 



3nR ± V48 + 1057r 2 i? 2 
Wher6 ' m± = 8(4^^ + 1) • (3 ' 5) 

As R becomes large, m± becomes smaller as 1/R and M* accesses all possible values ranging 
from small to large. Thus the nontrivial solutions of (3 g = and /?m 2 — gives dense lines of 
fixed points occupying a region. Depending on all these values of M*, g* ^ will have different 
values from ( |3.1| ). Note that all these nontrivial fixed points corresponds to c = 1 as they have 
the string susceptibility exponent of c = 1 (70 = 0) irrespective of the different values of g* and 
M*. It can be mentioned here that as R — > 00 the magnitude of g* becomes smaller and smaller 
and eventually all the pair of fixed points for ±g* 7^ coincides with the Gaussian fixed point. 

Now let us look at the general shape of the trajectories flowing to these c = 1 nontrivial 
fixed points g* 7^ 0. The RG equations /3 g — and (3m 2 — can be combined to give 

8M 1= 2 !)i F 32 -MHF^F 32) ) 
dg -1 + 9 2(2F 93 - 3F„) 

For large R this equation takes the following form where the variables are easily separable 

dM (F g2 -M 2 (F g2 + F g2 )) 



dg gM(2F g3 - 3F g2 ) 

Using the explicit functional form at large R this becomes 

dM M(8ir 2 R 2 M 2 - AttRM - 8M 2 + 9) 
~dg~ ^(247r 2 i? 2 M 2 - IOttRM + 3) 



(3.7) 



(3i 
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The leading behavior (considering 1 /R as small parameter) indicates, 

M 2 -/ 6 . (3.9) 



dM dg _ , r2 



In the last section we will come back to this behavior in the context of Neumann and/or Dirichlet 
boundaries of 2D string theory inserted by the integration of the fundamental fields (quarks) 
v* and v. For Dirichlet boundaries one expects large value of M*. Then the boundary fields 
become uncorrelated in time. One can consider the scaling function of M 2 as the boundary 
cosmo logical constant fiB- Since Hb is expected to scale with the bulk cosmological constant \i 
as 

/ifi^V^) ( 3 - 10 ) 

and since the scaling function of g gives the renormalized bulk cosmological constant, the above 
scaling behavior in (|3.9|) is nothing but an indication of the presence of various boundaries of 
the 2D string theory. 

3.1 The critical exponents 

Let us now go back to the matrix partition function we started with. After completing the RG 
transformations, it obeys the relation 

Z N+ i[g,M,R] ~ [\(N,g,M,R)] N2 Z N [g' = g + 5g,M' = M + 5M, R = R + 6R] . (3.11) 



This leads to the Callan-Symanzik equation 
for the string partition function (or the world-sheet free energy) 



F[g,M,R]^r[g,M,R]. (3.12) 



T[g, M,R} = ^ In Z[g, M, R] , (3.13) 

with 

7 = 2. (3.14) 

The singular part of the world-sheet free energy T s is given by the solution of the homogeneous 
Callan-Symanzik equation. The inhomogeneous part defined by the change in the prefactor A, 
contributes to subtleties in the free energy. 

Let us now discuss the critical exponents for the scaling variables, the renormalized bulk 
cosmological constant A = 1 — g/g*, the renormalized mass (or in the context of boundaries, 
the renormalized boundary cosmological constant) M. = 1 — M/M*. Introducing the matrix 
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= dAi ' ^ ^ 

the homogeneous part of the Callan-Symanzik equation, satisfied by the singular part of the 
free energy, can by rewritten as: 



d d d d 



F a [A,M,R]=0, (3.16) 



where, Q{S are eigenvalues of the matrix flk,i- They are nothing but the scaling dimensions of 
the relevant operators. The general expressions of Qk,i s an d the eigenvalues Q l5 Q 2 , for different 
fixed points (both Gaussian and the nontrivial one) for a general nonzero h are evaluated in 
the Appendix B . For h = 0, the (3r term drops out from the Callan Symanzik equation. The 
singular behavior with respect to the renormalized cosmological constant goes as, 

T s ~ A 2/f7l /i[^ ni A ] f 2 [N Q ' 2 M ] . (3.17) 

Comparing the above expression of T s with the matrix model result T s ~ A( 2 ~ 7 °) f[N 2 ^A}, 
or using the standard definition of the susceptibility T ~ f^l m=o ~ A~ 70 , the string 
susceptibility exponent 7 is given by 



7o ~(2-2/fti). (3.18) 

Note that in our analysis 2/71 ~ Qi, i.e. 71 ~ 2/Qi is consistent with the matrix model relation 
7o + 7i — 2. This relation is independent of the explicit values of 70 and 71 and is easily 
obtainable from the consideration of the torus. The string susceptibility exponent at genus G 
is defined by 

1g = 70 + ^71. (3.19) 

Referring to the Appendix B, we observe that, for our nontrivial fixed point, ill = 1 — 5h — 1 
(h — 0) , and hence 70 = 0. This shows that our nontrivial fixed point is a c = 1 fixed point 
at any R for h = 0. This c = 1 nontrivial fixed points are repulsive with respect to the flow 
of the parameter g while the gaussian fixed point is attractive. Similarly one can analyze the 
critical exponents for the scaling function Ai, the boundary cosmological constant, which will 
characterize different boundary conditions. Such an analysis will be useful to understand all 
possible boundaries and to identify the Neumann and Dirichlet case. 

Note that in this RG, typically the nontrivial fixed point is always situated close to the 
Gaussian fixed point. Hence, flu, ^22, ^21 components of the scaling dimension matrix are 
small. Thus Qn ~ Qi and Q 2 is also small. 
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3.2 The logarithmic scaling violation of the singlet free energy 

In the previous section, the behavior T s = A 2 of the singular part of the free energy near the 
c = 1 nontrivial fixed point (70 = 0) as a function of the renormalized couplings or the scaling 
variables is consistent with the continuum prediction, the KPZ-DDK scaling law |5], |6], [7|. This 
power law dependence on A is present in all known c < 1 theories. However, for c = 1, matrix 



model predicts a logarithmic deviation to the usual power law scaling [^, [23|, |24], |25], . This 
is known as the logarithmic scaling violation of the c = 1 matrix model. In terms of the singlet 
free energy or the ground state energy E , 

~N^ S = E = -iW-^- + . . . , 

In A 

N 2 1 
or, E = —- — \i In // + ..., where, A = — — /iln/i. (3.20) 

<±7l ' 271 

In the inverse Laplace Transform with respect to the area A, the sum over surfaces of fixed area 
A behaves as 

The continuum methods seem to be insensitive to this peculiar dependence. 

These scaling violations of the c = 1 matrix model can be naturally explained considering 
the lack of the translational invariance in the Liouville dimension of the 2D string theory, 
that is manifested by the <fi dependence of the background fields and is needed to maintain 
the conformal invariance. Here, the Liouville dimension arises indirectly from the semiclassical 
dynamics of the matrix eigenvalues while the Euclidean time of the matrix quantum mechanics 
provides the other dimension. The logarithmic factors in the expression of energy actually 
reflects the logarithmic divergence of the Liouville volume as the critical point is approached. 
The scaling violation was argued to arise from the unusual dependence of the tachyon potential 
on the Liouville field, T(0) ~ p7 |. Later it was shown that [[5^, considering interactions 



that represent touching random surfaces give rise to new critical behavior. This corresponds 
to the other branch where tachyon potential has the ordinary Liouville form, T(0) ~ e^, with 
simpler scaling 

J- s (A)~A 2 lnA, JT(A)~J_. (3.22) 

We would now like to see how one can get this behavior peculiar to c = 1 from the Callan- 
Symanzik equation. The usual scaling violation, T s ~ A 2 / In A, directly follows from the 
solution of Callan-Symanzik equation ( |3.12| ) near the c = 1 nontrivial fixed point 



\r 9 . d 

N A — 

dN dA 

up to a term that vanishes in the scaling limit. 



iV 2 ^ s = 0, 
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Also it is observed that the renormalization of A giving rise to the right hand side of the 



inhomogeneous Callan-Symanzik equation (|3.12|) can contribute logarithmic terms in the be- 
havior of the singlet free energy [[37]. The reason right hand side is capable of contributing to 
the free energy (universal physics) is that it actually measures the change of the world sheet 
free energy with the change of scale and can be expressed in the same footing as the terms in 
the left hand side of the Callan-Symanzik equation 

dT dT 



dX 



(3.23) 



Thus, using our Callan-Symanzik equation ( |3.12| ) and the linearized beta function equations 
( |2.35|) around the nontrivial fixed point ( |3.1|) we have , 



,r 9 A d 

A— 

dN dA 



N 2 T = 13, 



d(N 2 T) 
' dX 



-AF g3 (R)g* + A 2 ^f 



This gives 



2T- A 



dT 



dA 
d{T/A 2 ) 



N 



-2 



N 



-2 



AF g3 {R)g* + A 



F g3 (R) F g3 (R)g* 



2 F g3 



2A 



A 2 



dA 



(3.24) 



(3.25) 



Integrating both sides gives rise to the scaling 



E = —N 2 T 



-~F g3 (R) A 2 In A. 



(3.26) 



However, this dependence does not give rise to the logarithm in the inverse Laplace transform 
with respect to the area. This is rather the behavior corresponding to the other branch of 
the tachyon of the Liouville theory. Perhaps this reflects the fact that even though we are 
not explicitly considering the branching interactions representing the touching random surfaces 
generated by the redundant higher order terms, they have indirect effect in the change of the 
world sheet free energy. 

Note that in this analysis we are not looking at the M. dependence of the free energy since 
typically Q2 is very small and hence can be ignored in ( |3.12|) or in ( 3.16|) . Alternatively, one 
can also consider a fixed mass parameter M in doing the analysis. 



3.3 The Radius Dependence of the free energy 

As we have discussed before in context of the analysis of the structure of the beta function 
equations, one can observe an indirect impact of the world sheet vortices on the free energy in 
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the sense that it flips its sign at the self-dual radius R = 1. This radius is known to be the 
inverse temperature for the Kosterlitz-Thouless phase transition due to the liberation of the 
world-sheet vortices . To observe this we would demand that the mass parameter M is fixed 
to a certain value (for simplicity we choose M — 1) while (3m 2 is set to zero. Then for nontrivial 
fixed points of (3 g = 0, the parameter h can be determined as a function of R by solving the two 
simultaneous equations (3 g = and Pm 2 \m=i = for the two unknowns h(R) and g* . This 
gives, 

h(R) = - v . (3.27) 

(4 F g3 (R)/F g2 (R) - l) 

In low R approximation h(R) ~ 1/R 2 and thus the change of the world-sheet free energy around 
the nontrivial fixed point becomes (considering g* to be very small) 

(3^\{l/R 2 -l). (3.28) 
This indicates a phase transition at R = 1. 

h-1 
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Figure 3: Behavior of h(R) — 1 at relatively large R. 



However, a plot of f3\ « \{h{R) — 1) using a simplified h(R) ( h(R) ~ ^T+TW/r ) m 
relatively large R expansion shows the discontinuity at R = 1.2 (Figure^). In this discussion the 
logarithmic constant in the equation for (3\ ( |2.35| ) is ignored as it merely indicates some vacuum 
normalization and can be conveniently absorbed. Another important point to be mentioned is 
that the high R approximation for h(R) is not quite efficient to compute the correct exponent 
of the c = 1 fixed point (in the sense that h(R) does not go to zero for large R instead approaches 
— 1 in that limit). Since this is purely dependent on the contribution of the single line Feynman 
diagrams in integrating out the quark loops, we need further improvement in the computation 
of the diagrams. 
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4 The T-duality of the singlet partition function 



In this section we would analyze that in what sense our RG can capture the T-duality respected 
by the singlet sector of the MQM. Let us consider the matrix partition function in the eigenvalue 
representation by diagonalizing the matrices by SU(N) transformation and suitably integrating 
out the angular degrees of freedom. The diagonalization of $(£) gives, 

N 

M*) = ■ (4.1) 

k=l 

Hence, 

N 

Tr $ 2 = + E( A * - A i) 2 l4yl 2 > = (n f (W))* • ( 4 - 2 ) 

Here the gauge field A acts as a lagrange multiplier that projects onto SU (N) singlet wave 
functions which depend on the N matrix eigenvalues only. In the recent works on N x N 
matrix quantum mechanics as the quantum mechanics of string tachyon field on the N unstable 
DO-branes, this gauge field is identified as the non-dynamical gauge field in the open string 
spectrum corresponding to the vertex operator t. Due to the periodic boundary condition on 
the $(£) and on the Q(t)-s the gauge field Ay(t)-s are constrained: 

fexpi / A(t) dt = I. (4.3) 
Jo 

Taking into account of the above constraint (which actually projects to the singlet sector) in 
the measure of the path integral, the (singlet) partition function can be formally written as (see 
the review JJ): 



Z N (g,R) = 4°> / I]l'Mi^ r *" i:f ' 1 ' ; ™ , ^ & "^ 1 , (4.4) 

r J k 

where, the cIr and are the dimensions of the R-th representation of U (N) and that of the 

(2) 

subspace of vectors with the zero weights in R and C r is the value of the quadratic Casimir 
operator in the r-th representation. 

CP ~ N . (4.5) 

Note that, in this process, all the Vandermonde determinants are gone but the fermionic 
statistics is still maintained for the eigenvalues are now periodic: 
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A 1 (0) = A j (2tt j R), i = l,...,N. (4.6) 

In a similar sense to our scheme, let us now begin with the partition function of (N + 1) 
eigenvalues and observe the effect of integrating out one of the eigenvalues in presence of the 
rest of the N eigenvalues. Physically this would mean integrating out one of the fermions in 
the presence of the interaction due to the rest of the N fermions. Then the Fermi sea readjusts 
it's height and this renormalizes the parameters of the theory like the cosmological constant. 
In recent works a similar picture has been discussed in context of closed string radiation due 
to the rolling open string tachyon on unstable DO branes represented by the unstable matrix 
eigenvalues. We will come to that in the next chapter. 

Now the partition function of (N + 1) eigenvalues can be split as, 

r J k=l 



/ 



Let us consider the simplest form of real (Hermitian Matrix model) and periodic eigenvalue, 

Ajv+i(0 = Ajv+i cos t/R . (4.8) 
Then the measure over the (N + l)-th eigenvalue trajectory (to be integrated out) becomes, 

N N N 

V\ N+ i(t) = l[d\ N+1 cos(2irk/N) = l[dX N+1 (l-0(l/N 2 )) ~ ]JdX N+1 (4.9) 

k=l k=l k=l 

Hence the integration over Ajv + i gives, 

N 



In+i = 1 I / d\ N+1 e 



1 X 2 

2 JV+1 



lt R dt ( sin2 t/R+cos 2 t/R) ^ R dt cos 3 t/R 



n/ 

fe=i J 

1 [y^JV 1 rl-nR dt i 

e 2(AT+1)2 L^j = l (\ N+l -\ j ) 2 JO cos2t/ijJ 

N 

" II / dXN +i e -^ N+1 X R+1 / R ^ = (2/(N + 1)(R + l/i?)) f . (4.10) 
fc=i •* 

As before, considering the rescaling of the fields 

Ai(*) - pAj(f) , i = l,..,N, (4.11) 

and setting the coefficient of \X^{t) to one (which in this simplest case simultaneously sets the 
coefficient of |A^ 2 (t) to one too, and the mass term does not run with the scale change) we have, 
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P=(l + ^)^-(l-7^ + 0(l/iV 2 )) 



N' v 2N 
Hence, the partition function can be rewritten as 



(4.12) 



w+i 



(g,R) 



1 



. (N + 1)(R + 1/ R) 

r N 

J^drd^ j\{V\' k {t)e 

1 , 



i^j (X'.-X'.j 1 



k=l 



9 = 9 



9 
2N 



N 2 ' 



C = 1- — + 0(At' 
N K N 2 ' 



(4.13) 



This leads to, 



'N+l 



(g,R) 



z N (g',R) 



(N + 1)(R + l/R) 
1 / 2 



(i-^ + o(Viv 2 )) 



N 



(1 - l/NY 



n N 2 



(4.14) 



-ATW \R+1/Ry 

Note that, due to the cosine, which is the simplest choice for the eigenvalue Ajv+i, both the 
interaction terms corresponding to the cubic self- interaction of A n + i and the mutual (repulsive 
coulomb) interaction with the rest of the N eigenvalues drop out in the integration In+i- As 
a result in the equation (|4.7| ), the integration over the + 1-th eigenvalue becomes just a 
overall prefactor (contributing to the change of the world-sheet free energy associated with the 
readjustment of the Fermi level in response to loosing one of the fermions) as in ( |4.13| ), unlike 
the more general situations studied in the previous chapters where the integration over a part 
of the matrix degrees of freedom does produce nontrivial one loop correction terms that adds 
up to the left over A^-fermion partition Z^{g' ', R) besides producing the overall prefactor. Even 
in simple eigenvalue representation, a more general choice of the functional form of Xn+i would 
produce such correction terms. In those cases the beta function for g, (3„ 



*(¥) 



will have higher order terms in g (which can give nontrivial fixed points). Note that the leading 
term in the beta function, — |, is the same in any case. 

Now the equation ( |4.14j ) implies the Callan-Symanzik equation for the world sheet free energy 
r\g,R] = -falnZ\g,R]: 



N m-^g +1 . 



T[g,R] 



X i f 2 



)+0(l/N), 7 = 2. (4.15) 



In terms of the scaling variable A 
scaling limit is given by, 



g/g*), the Callan-Symanzik equation in the Double 



\r 9 1 A 9 

N A h 7 

ON 2 dA \ 



JF[A, R] ~ -ln( 2 -^). 

1 ' J 2 y R+l/R J 



(4.16) 
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Solution of the homogeneous equation is same as before. Using the same ansatz as before, 
JF[A, R] = N- 2 f 1 (A)f 2 (R), and solving the inhomo geneous part we observe that the world 
sheet free energy has a functional form of f{R+ i/R) with respect to R: 

F[A,R] = -]nA ln( ^ + 2 ). (4.17) 

This shows that the large N RG is capable of capturing the T-duality property respected by 
the singlet partition function. However, two comments are note worthy here. First of all, in 
this simplest case, the beta function does not receive higher order corrections in the coupling 
constant. Hence the fixed point is a trivial Gaussian fixed point with inadequate critical expo- 
nent (fii = — ~,7o = 2 — 2/fix = 6) to describe the c = 1 fixed point (70 = 0). Actually, from 
the results of the more general set up in the previous chapters, this Gaussian fixed point could 
be thought of to be overlapping with a pair of double zero of the beta function at infinitesimal 
distance from it for large R (where the singlet free fermion picture is meaningful). With a more 
general choice of the functional form of the eigenvalue, the gaussian fixed point will be resolved 
into a nontrivial double zero of the beta function with appropriate critical index to describe 
c = 1, as obtained in the previous chapters. This c = 1 fixed point would explicitly show T- 
duality in exactly the same manner from the gaussian part of the integration over the N + 1-th 
eigenvalue in presence of the rest of the N eigenvalues. Secondly, we note that the T-duality 
arises in the usual sense from the kinetic and the quadratic term of the action in the wave func- 
tion of the N + 1-th fermion. This wave function can be compared with the free fermion wave 
functions for the singlet sector of c = 1 on large circle. A small observation is that in the more 
general set up (with M = 1 as in here), the gaussian part of the integration over a part of the ma- 
trix degrees of freedom produces the functional form /(sinh 7iR/ii 2 R 2 ) ~ /(^ + ^ + 0(i? 2 )) (to 
compare see equations fl2.34|) ), instead of producing f(R+l/R). This shows how the T-duality 
is broken at this level. 



5 Comments on D-brane Decay and Rolling Tachyon 



The recent observations 3^, |4"T 
dynamics of D-branes in c 



J show that the matrix model can be realized as the effective 
1 non-critical string theory. We briefly review the basic picture 
before going into the matrix model RG interpretation of the open string rolling tachyon in string 
theory. 



5.1 Review of D-brane Decay in 2D String Theory 

According to this conjecture, the matrix $at itself can be seen as an open string tachyon field 
on .DO-branes and the matrix potential as the tachyon potential. The SU(N) symmetric matrix 
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quantum mechanics in an upside down inverted harmonic oscillator potential, i. e. in the double 
scaling limit, exactly describes open strings on N DO-branes and is dual to Liouville theory 
coupled to c = 1 matter, that describes the two dimensional closed string theory together with 
its DO-branes. The .DO-branes, on which the tachyon field resides, are described by the localized 



boundary states for the Liouville coordinate <p introduced by A.B. and Al.B. Zamolodchikov |45 
tensored with Sen's rolling tachyon boundary states with Neumann boundary condition for the 
free time direction t . They are localized in the strong coupling region <fi —>■ oo far from the 
bulk region of large and negative <fi (signalled by the absence of "bulk poles" in the one-point 
function of the closed string vertex operator on the disc). These branes are parametrized by the 
pair of integers (m, n). Among them only the set (1, n) has smooth behavior with no singularity 
in the classical limit h — > 0, where b measures the "rigidity" of the 2D surface to quantum 
fluctuations of the metric. The one labelled as (1, 1) contains operator that matches with the 
massless tachyon in the open string spectrum and is the only one which has been found to 
be consistent with the standard loop perturbation theory and consequently is identified as the 
boundary state at one-dimensional infinity or the "absolute" of Euclidean AdS2 {i.e. classical 



Lobachevskiy plane) |f45| . It is still not clear what the other boundary states with m^l, n^l 
correspond to. 



There is also another class of Liouville boundary states |55j, extended in the bulk weak 
coupling region of — > — oo (exhibited by the poles in the bulk one point function). They 
correspond to D-strings when the time direction is taken to be Neumann. The continuous 
families of such boundary states are parametrized by the uniformization parameter s, given by 

l^B = ^/yU COsh 7TS . (5.1) 

Here /is is the renormalized boundary cosmo logical constant and /i is the renormalized bulk 
cosmological constant in the Liouville theory: 



Hb = {-K^Bol{b 2 )) , 

7« = fgL, (5.2) 

in the sense that, in order to get finite amplitude, they are kept finite in the limit b — > 1, i.e. 
cl — > 25. All such D-branes in this family have a continuous spectrum of open string states, 
namely, the massless open string tachyons. For the choice of s = |(2n + 1), m G Z, one can set 
[ib = 0. The work of || uses the minimal values s = ±|. Increasing m gives rise to increasingly 
tachyonic open string modes. 
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As the tachyon is rolling down the potential, one can identify the boundary cosmological 
constant /is of the Liouville theory with the time dependent tachyon potential A cosht as follows: 

^s{t) = y/Ji cosh s(t) = A cosh t. (5.3) 

As t — > — oo, Hb — > 0, which corresponds to unstable D-string picture. As t — > oo, hb and 
s grow to large values. In this process, the unstable D-string one starts with, decays into 
closed string vacuum of Liouville theory with its localized DO-branes and all the open string 
excitations on the D-string are pushed to (f) — > — oo. In the matrix model picture, an unstable 
brane corresponds to a free fermion on the top of the inverted harmonic oscillator potential 
rolling down as an unstable eigenvalue z{t) = Ae* from the top to the Fermi level 

The amplitude for such D-brane tachyonic decay with momentum P and energy u = \P\ is 
given by the expectation value of a normalizable bulk operator t^p = exp[(2 + iP)cf) + i\P\t] 
inside the disk. As was shown in [[KJ this disk one-point function agrees with the outgoing 
radiation amplitude derived from c = 1 matrix model (taking into account the leg-pole factors). 
The disk one-point function in the two dimensional string theory is given as follows: 



A(u,P) = A t {u)A L {P) = (v u \B x )(v P \B s ) , (5.4) 

where the two factors denote the time part (obtained in and the Liouville part (obtained 
in |55[ and [45] for the two kinds of branes) respectively. Although the end result of tachyonic 



decay for both the branes are the same and the time part of the amplitude is identical, the 
Liouville part or the bulk one point functions for the two cases differs with the absence of pole 
for the localized brane. 



The time part of the amplitude A t {ui) = {v u \B>^ depends on the choice of contour [56[ and 



is described by the boundary interaction^] A cost/a'. Note that, in our convention, t denotes 
the Euclidean time. In the Lorentzian coordinate X° = it, the boundary interaction is given by 
A cosh X°/a'. The choice of Hartle-Hawking contour gives rise to the following result 



. , . . r . ,. 7rexp[— iy/o/u log X] _ \ 

A t (u) = (expM) = P 1 v & J , (5.5) 

sinh 7tv a co 

having pole at u — 0, signaling the presence of continuous spectrum in the on-shell open string 
channel |57|, H- 

The Liouville part of the amplitude for Dl-brane decay (for the choice of s = ±z/2) can 
be evaluated from the on-shell CFT calculation of the boundary state constructed by Fateev, 
Zamolodchikov and Zamolodchikov [|5f| resulting the total amplitude to be 



1 Here A = sin7rA in order to match with the trajectory considered in [Em. 
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A D1 (u,P) = 2n 2 e~ iS ^ e -^ p ^ X costtV^P (5 g) 

sinh 7tv a'o; sinh7rv ck'P 



The corresponding phase factor is given by 



T(-iVa'P) 

The corresponding disk one-point function for the DO-brane has no pole [^UJ. The Liouville 



part is described by the Zamolodchikov boundary state ESI and the total amplitude is given by 



A D0 (uj,P) = -sgn(P)7r e- iS{p) e -^W°s\ (5.8) 

Let us now try to understand this from matrix quantum mechanics picture. The physics of 
the matrix model in eigenvalue representation is that of iV non-interacting fermions, moving in 
the potential V(X). In the double scaling limit, the potential is that of an inverted harmonic 
oscillator. Each fermion occupies a volume of 2n/(3 = 27rg 2 /N in phase space, where 1/(3 plays 
the role of h. All levels with single-particle energy |t> 2 + V < /xp are filled. As g increases or (3 
decreases, the phase space volume occupied by the fermionic states grow. At the critical value, 
g — > g c , the fermi level reaches the top of the potential, /xf — > /i c , where the fermions start 
to spill over the barrier. The free energy is singular at this point. In the double scaling limit, 
iV — > 00 and g — > g c , holding the string coupling g s = ft^ 1 = [/3(/j c — /ip)] -1 fixedQ. In other 
words, the ratio of the gap between fermionic levels (which is of order /?, ~ 1//3) to /i = /i c — fip 
remains constant in the double scaling limit. In this limit the physics is that of the Fermi sea 
in inverted harmonic oscillator with the hamiltonian H(v , A) = (v 2 — A 2 )/2 + /i, where A, v are 
canonical conjugate variables. The collective motions of the fermions are classically described 
in terms of a time dependent Fermi surface that separates the filled and the empty regions in 
the phase space. The fermions on the surface moving freely in the inverted harmonic oscillator 
potential are described by, 

D t v = \, D t \ = v, (5.9) 

where D t denotes the co-moving derivative following a phase point on the surface. Hence the free 
fermions on the Fermi surface execute simple hyperbolic orbits. However the time dependent 
profile of the Fermi surface is more complicated. For small perturbations, the profile can be 
described by the positions v±(X, t) of the upper and lower surfaces of the Fermi sea at each time, 
satisfying, 

d t v±(\,t) = \-v±(\,t)d x v±(\,t). (5.10) 

2 In terms of A = g - g c ~ fi \nfi, we have kept g s = [TVA^ 70 )/ 2 ]- 1 = (NA 1 ^ 11 )^ 1 fixed in our RG analysis 
in previous sections. 
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Therefore, the static Fermi sea satisfying the equation of motion ( 5.1 0|) is given by, 

v±{\) = ± v / (A 2 -2/i), i.e. ^(v 2 - X 2 )+p = 0, (5.11) 

which also given by the ground state of the hamiltonian H(v,X) = (v 2 — A 2 )/2 + p . For 
simplicity, we assume that, in the ground state, the local maximum of the potential is at zero 
energy level (/i c = 0), and the energy at the Fermi level is fip = — /x, measured from the local 
maximum. Identifying 



v(X) = — = v/(A7«' - 2fi) (5.12) 

as the velocity of the classical trajectory of a particle at the Fermi level, we define a new spatial 
coordinate r, the classical time of motion at the Fermi level. 

A single DO-brane state of the continuum theory corresponds to a single eigenvalue or a 
fermion excited from the Fermi surface to higher energy, e.g. to the top of the upside down 
potential. This corresponds to putting the Neumann boundary condition on the world sheet 
field X(ai,a 2 ). Turning on an exactly marginal boundary operator that interpolates between 
the Neumann and Dirichlet boundary conditions on X{a%, 02) makes the unstable brane to roll 
down to D-objects localized at large Liouville direction <ft. This corresponds to considering an 
unstable eigenvalue as the probe fermion decaying in the presence of the N free fermions, as it 
executes classical Euclidean motion in the forbidden region along the trajectory 



A(t) = sfioip X cos(t/a') . (5.13) 



In order to relate to the picture of the rolling tachyon in fl£6| , the rolling boundary state 
parametrized by A, where sin7rA = A, is identified with the decaying eigenvalue starting at 

A = -y/2oFp, A, (5.14) 
where the energy of the state is given by 

E = p cos 2 7T A . (5.15) 

The time delay in the classical evolution of the trajectory relative to the classical trajectory at 
the Fermi level is therefore 

f T dr' f dX , 
-At = - —= = / —. ~ In A , (5.16) 

J y/a' J VA 2 -1 

for large A. The parameter A G [0, 1]. The case A = 0, i.e. A = 1, describes an eigenvalue at 
the top of the inverted harmonic oscillator potential (V = 0, E = —p) and X — 1, i.e. X = |, 
describes and eigenvalue on the Fermi sea (V = —p, E = 0) . 
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To relate such physical process in the free fermion theory to that of the D — 1 noncritical 
string theory in the dilaton and tachyon background, one needs to bosonize the non-relativistic 
free fermions in a suitable way. One way to view this is to express the small fluctuations of the 
time dependent Fermi surface in terms of the massless scalar X(t, t), where r = In A ( |5.16| ) fT3 |. 
In this respect equation (|5.10|) can be generalized to the time dependent case as 



where 



v± (X,t) = T A±-e ± (r,t) 



e ± (r,t) = ±d t X(r,t)-d T X(T,t). 



(5.17) 



(5.18) 



The right and left moving fluctuations e± (r, t) are proportional to the right and left moving 
currents, J R =: ip R ipR : and Jl ='■ iPl^l : > generated by the chiral fermionic variables as given 
by the bosonization relation 



'Ik 



2n 



: exp 



exp 



i^K I dr\d t X - d T X) 
dT'(d t X + d T X) 



W7T 



(5.19) 



(This is because {Jl + Jr) = d T X and (J L — J R ) = ^ d t X) . The second quantized 
hamiltonian for the free fermions is written in collective fields t) = ^ aiipi{\)e~ teit , where 
ipi are the single particle wave function and are the corresponding annihilation operators. One 
can expand the fermionic field ip(\,t) in ipR and ipL as 



A p\ 

cx])[-i I dX'v(X') + iir/i) 4> L (\,t) + exp[i / dX'v(X') - iir/A) ip R {X,t) 



and rewrite the second quantized hamiltonian 



(5.20) 



H 



l,v" 



dr 



i(ip R d T ip R - ip{d T ij L ) + —(d T i)ld T il) L + d T ip R d T ip R ) 



5v 



f'2 



A y v 3 2v i 



(5.21) 



Here v' = dv/dr. Solving the equation ( p,16| ) for A , A = cosh(r/va') ~ expr, r — >• oo . 
Hence v = i/2/Isinhr ~ -^/S/Iexpr. This shows that the 0(l/v 2 ) terms are not negligible in 
the double scaling limit where we keep fx fixed, unless we go to the asymptotic or the r — ► oo 
region. This maps the second quantized hamiltonian of iV non-relativistic free fermions to the 
two dimensional Dirac hamiltonian. This shows that r is the natural spatial coordinate in 
terms of which the fermionic system has a standard Dirac action to leading order in (3 . Also in 
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terms of the bosonization variables the hamiltonian in the r — > oo limit takes the form of the 
canonically normalized free scalar hamiltonian 



H~~jdr [{d t Xf + {d T Xf + e~ T 0(X 3 ))] , (5.22 



with the equation of motion (df — d^)X = e~ T 0(X 2 ). This also shows that r can be identified 
with the zero mode of the Liouville field . One can restrict A to lie between the two turning 
points of the classical motion ( ±y/2a'p, ) , i.e. A G [—1,1] , or equivalently restricting r to 
lie between and T/2 where T is the period of classical motion {§ dr = 2i\T\ One then 
imposes a standard bag-like boundary conditions on the chiral fields to get a consistent picture. 
Introducing the cut-off on r, the energy levels are approximated by 

e„~ -fi + n/T + 8'(n/T, /i) + ... (5.23) 

where the last two terms can be dropped in the limit T — > oo. 

Thus, if a single fermion starts near the top of the upside down potential it will roll down 
to infinity. In the region r 3> 1 , A ^> 1 , it becomes approximately relativistic (since u> 1). 
The final state could be expressed in terms of the chiral fields bosonized as 

fa, R {T>t)\ty = e i * / * XLMT ' t) , (5-24) 
(where, Xi + Xr = X ) . Since the bosonization scalars have the mode expansion 

J 27tV Ah L J 
a final state of the decay would be given by 

J dr^(T)ij L \0) = J dT^iry V5F*i(T, *)| ) = J dr^^e 2 ^ 1 A e ~ ipr a [ p \0) . (5.26) 

Thus keeping in mind the form of the coherent state obtained in the tree level string theory, 
\ip) ~ exp ( J a) p A |0)), one can extract the decay amplitude from the bosonization scalars 
as 

A ~ i2\fKe~ iEr , \p\=E. (5.27) 
This is the same as the amplitudes calculated from the continuum method. 

5.2 The RG Picture of Tachyonic Decay 

One can view such a process of decay of an unstable eigenvalue describing the open string 
tachyon roll-down to 2D closed string vacuum, as a Large iV Wilsonian RG in the singlet sector 
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of the matrix quantum mechanics leading to the c = 1 fixed point (in a similar fashion to our 
calculations in the previous sections). The RG equations are nothing but a statement of response 
of the free fermion system as successively one of the fermion is integrated out in the presence 
of the rest of the free fermions. Hence, following our calculation, by performing the Large 
N world sheet RG for the singlet Matrix quantum mechanics in the eigenvalue representation 
(with Ajv+i(£) = z(t) = yjla'ji A cost/v^o 7 ), one can compare the change of the world sheet 
free energy in the Large N RG with that due to the closed string emission picture. One can 
do this comparison both for the singlet MQM on infinite line and that on a large circle. In 
the later case one would have to consider closed string emission amplitude for compact time on 
large circle. This comparison, on one hand gives rise to a physical picture of the Large N RG 
we considered in this paper, on the other hand it provides a way to test our results. Note that 
the desired c = 1 end point, corresponding to the 2D closed string vacuum, always emerges in 
a very nice and robust way in the RG calculation with all it's critical properties. 

We would now try to understand the emergence of this amplitude from the point of view of 
the world sheet RG by turning on an unstable eigenvalue 

z(t) = y / 2a'fL A cos(t/a) , (5.28) 

and integrating it out in the presence of the rest of the N free fermions. The Fermi sea, in 
response to this, would readjust its height infinitesimally which is reflected in the renormalization 
of the couplings. Later we will focus on the more general context of the Wilsonian world sheet RG 
of the (N x N) arbitrary matrices on circle rather than working in the eigenvalue representation 
only, which is good for the large radius. 

The partition functions before and after turning on this unstable probe eigenvalue respec- 
tively, are given by, 



N 



Z N [g,t]= / HdX^t) Ajv(t) exp 
J i=i 

N 

Z N +i[g,t} = / Y[d\(t) A N (t) exp 
J z=i 

J dz(t) exp |^ 



NT* 1 • 

-- 2 J dt'i-Kttr+viW))] 



9 Z 



dt'[-\ t (t') 2 + V(\ t (t'))} 



( Af + I) ft 1 i 1 „ 

' < dt'{-z 2 (t') + V(z(t'))}\ exp [-#*(*))] , 



(5.29) 



where the operator 



1 



- (t>{z{t)) = log 



Ajv + i(t) 
A N {t) 



N 



^ log Xi(t)-z(t) , N~l/h. 



(5.30) 



i=l 
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The Laplace transform of the operator (j>(z(t)) can be seen as the loop operator 0(1) that creates 
a closed (Dirichlet) boundary of length I on the world sheet. Its expectation value gives the loop 
amplitude [|58|| . In the limit of the vanishing boundary length, the operator becomes tachyonic 
and its excitations correspond to the excitations of the closed string tachyon: 



/oo /*oo 
dqe~ iqt dl e~ lx{T) 0(l,q), 
-oo JO 

0(1, q)~~f dt ^ J dT e ~ lX(T) d - X ■ 



(5.31; 



The mode expansion of the massless scalar field X is related to the amplitude A of the closed 
string decay in the sense of final state of the evolution being somewhat like a coherent state 
hj, as shown in ( |5.26|) . Thus in the matrix model, one can also compute the decay amplitude 
by computing <p(r,t). (The reference |4l| describes a similar thought to relate 0(r, t) to the 
amplitude. However, they chose <p(r,t) to be in the same footing as X, instead of fl5.31| )). 



Let us now consider the partition function for the N + 1 eigenvalues (|5.29|) and rewrite the 
right hand side in terms of the renormalized coupling g' 



^^dt'ilut'f+viMt'))} 



nN 
d\i(t) A N (t) e 
i=i 



dz(t) e -^rdt'\\z\t') + v{z{t'))\ g i 



(5.32) 



where g' = g — ^ + 0(1/N 2 ). As before, extracting the partition function of iV eigenvalues, 
we have 



[9 A 



z N [g',t] 



dz(t) exp 



' I dt'{-z\t') + V(z(t'))}\ exp[-0(^))] 



if 



12 



1 



^(iV,^)(exp[-0(z(t))] 



where ipt(N,g') is given by 



MN,g')= I dz(t) exp[-^ J' dt'{h 2 (t') + V(z(f)} 



(5.33) 



(5.34) 



Then from the RG point of view the change of the world-sheet free energy (8T = log Z^+i [g, t]/N 2 - 
log Z]sf[g', t]/N 2 ) in the process of D-brane decay, given by the ratio 



Z N+1 [g,t] r 1 



Z N [g',t] 



MN,g') <exp [^<f>(r,t)]) 



(5.35) 



encodes the contribution from the amplitude A through a knowledge of </>(r, t) . In terms of 
the discrete Callan-Symanzik equation (|2.6|) , the left hand side indicates a flow from g — > g' as 
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(N + 1) — > N. The right hand side computes the contribution of the amplitude A to the change 
of the world sheet free energy as an 0(1/N) effect. 



(5.36) 



Here ip t is the single particle wave function corresponding to the probe eigenvalue. Note that 
for the singlet MQM on a large circle , the ipt thus calculated by the Large N RG is a function 
of (R + 1/R) which is a manifestation of the anticipated T-duality property. 

We would now calculate the right hand side of the relation ( |5.35| ) . Let us assume that the 
rest of the N eigenvalues belong to the static Fermi sea forming a closed string background 
given by 

A t = -yj2a!Ji A, i = l...JV. (5.37) 
Hence from ( 5.30|) and (|5.28|) and considering A ~ exp(r/Vo') for r 3> 1 and also N ~ 1/K , 



(r, t) = In [x/2a'fx(l - cos -=)] + r/Vrt . 



Thus, 



[ ^4 = A(N,g') (v^S(l 



a' 



cos—^=)) N /exp r/h^fa' 



(5.3J 



(5.39) 



ZwWit] Vtf 
Here the t dependent part could be absorbed in the overall prefactor of the ratio. As in Q5.35 ) 
effectively ~ r/y/a' contains the contribution from the decay amplitude to the change of the 
world sheet free energy. 

We will now go back to the relation of with the bosonization field X (and hence A) and 
verify that the expression of (f> we obtained in ( |5.38| ) does lead to the standard answer for A. 
Using ( |5.31| ) and performing the delta function integration over t , 

1 



(r, t) =7i dr' 



This implies 



{d T «i>) 



\/2a'jj, (exp t/vol' + exp r'/ v a') 
7i(d T ,X) 



(d T ,X(T,t)). 



\/2a'fi (exp t jya' + exp r'/ vV) 



Inverting and solving for X, 

X = — [ dr' \]2a'Ji (exp rj \fa' + exp t' / \fa') (d T '(f>) . 

7T J 



(5.40) 



(5.41) 



(5.42) 



Considering \d T i(p) = 5[t —t), we have X = ^ \/2pL ev? . Comparing with the relation of 
the mode expansion of X to the amplitude A we have, 



A ~ — a/ 2/i e^ 7 ~ e 

7T 



In A 



(5.43) 
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The time delay ( |5.16 )is thus given by \At\ ~ In A which is consistent. Thus we conclude that 
in our RG, integration of one eigenvalue in the presence of the others describes the decay of an 
unstable DO-brane with open string tachyon attached to it to the 2D closed string theory (with 
its .DO-branes) with an amplitude A given by 



d „, , d 



ON ' y ^dg 



1 , . 1 



N— -13(g)— + 7 F\g,t] = — \nA + — \nip t , A ~ e ln \ 7 = 2. (5.44) 



N N 2 



In our set up it is clear that there is a flow, g — > g' as (N + 1) — ► N, corresponding to the decay 
process. However, as we have discussed in the previous chapter, the eigenvalue representation is 
too simple to explicitly produce a produce a nontrivial fixed point of the flow (13(g) = —g/2) from 
the beta function equation. The flow in the eigenvalue representation only gives the gaussian 
fixed point corresponding to the inverted harmonic oscillator. However, still one can say that 
the flow is hitting the c = 1 fixed point situated at the infinitesimal distance from the Gaussian 
fixed point. Working in a more general set up of MQM of general N x N matrices, one can 
explicitly show from the beta function that the endpoint of the flow is indeed a pair of c = 1 
nontrivial fixed point ((3(g) = —g/2 + 3F(g,R) g 3 /2) dual to the 2D closed string theory the 
D-brane decay leads to. Also it is shown that for R — > 00 or even for large finite R (where the 
singlet free fermion picture is meaningful) the pair of c = 1 fixed points are drawn infinitesimally 
close to each other and to the Gaussian fixed point such that they are overlapping. Thus it is 
consistent to think that even in the simple eigenvalue representation the flow corresponding to 
the D-brane decay does end at the pair of overlapping c = 1 fixed point in the infinitesimally 
close neighborhood of the gaussian fixed point. This is derived in detail in the previous sections 
and is the key observation of this paper. 

We will now discuss the decay process in the context of RG analysis of the general N x N 
matrices 0jv(t). Consider the usual partition function Z N+1 of the matrix field <f>N+i(t) that is 
decomposed into (pnif) and the (1 x N) and (N x 1) row and column vectors v*(t) and v(t) 
respectively, 



Z N+1 [g,M] = J V N2 Mt) V N v*(t) V N v(t) 



cxp 



- N Tr / dt [^(D(f> N f + ^M 2 <P 2 N - + Dv*Dv + M 2 v*v - g v*(j) N v 



(5.45) 



As usual, the covariant derivative is defined in terms of the non-dynamical gauge field A in the 
open string spectrum corresponding to the vertex operator t, 



D(p N = <9, 



Wn 



[A, 



>N\ 



(5.46) 



36 



The gauge field A projects onto the SU(N) singlet wave functions by acting as a Lagrange 
multiplier. 

Now, one can see Z N+1 as the dual large N description of the 2D string theory in the 
presence of the space-time filling Dl-branes, the extended Liouville boundary states tensored 
with Neumann boundary state for the (Euclidean) time. One can similarly understand Z^+i for 
the (Euclidean) time taken on a (large) circle. Motivated from previous works [|7L BR BH], the 



authors of |^0| pointed out that integration over the vectors v*(t), v(t) inserts into the random 
surfaces dynamical boundaries (quark loops) wandering in the time direction. Although their 
model is essentially the same as ( |5.45| ), in the former case the couplings in the flavour part are 
introduced by hand rather than being determined by the original closed string action due to 
the parametrization (|2.7|) . 

As we have discussed in previous sections, precisely this is happening in our RG. The (N+l) x 
(N + l)-th matrix element gives the position of the wandering boundaries, which for simplicity, 
have been set to zero in our analysis . The sum of the 'single line ' one loop Feynman diagrams, 
generated by the integration over the vectors (v*,v) (see equations ( |2.21| ) and ( |2.24| )) , inserts 
in the path integral different operators coupled to the dynamical boundaries. The relevant 
operators coupled to these boundaries are of the form -^(Tr0 n — g n v*cf)^ n ~ 2S, v ). The massless open 
string tachyon would correspond to the operators J dt e iqt v*v(t). For a model without the kinetic 
term for (v*,v), integration over the vectors generates the macroscopic loops, the boundaries 
with the Dirichlet boundary condition on time. This happens when both the parameters M 2 and 
g acquire large values due to tuning. As pointed out in [f|D| , tuning the parameters M 2 and g 
together with N one should be able to get a scaling model with two independent parameters fiB 
and /i, the renormalized boundary and bulk cosmological constants, given by the relation (5.1). 



In our RG analysis we arrive at an analogous relation (|3.9|) with usual scaling between \xb and \x 



(with an accuracy of 80 %). This presumably indicate the presence of various boundaries in 2D 
string theory. However to realize whether these boundaries correspond to Neumann or Dirichlet 
case needs further investigation. Here, tuning of the matrix coupling constant g gives rise 
to a renormalized bulk cosmological constant /i. The scaling variable A = (1 — g/g*) depends 
nontrivially on \i 

A~— ^mAi, (5.47) 

Z7T 

as they are mutually coupled through the Fermi-Dirac distribution of the fermions in the grand 
canonical ensemble |p8| . On the other hand one can think of tuning the mass parameter M 2 
to control the boundary cosmological constant /i^ as 

fi B ~ 1 - M 2 /M* 2 . (5.48) 
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Hence solving the equation 



Pg . . dg 



A, SM "»,M,T), (5 ' 49 > 
the RG flows g = g(M, T) in the Double scaling limit near the nontrivial fixed points is expected 



to give rise to relations analogous to the trajectory (|5.1|). In equation ( |3.9|) we have solved this 
trajectory for large R (for matrix quantum mechanics on circle) and got a consistent scaling. 
We expect the fixed points of large M* obeying such scaling could be Dirichlet boundaries. 
However, since in our analysis we are dealing with bosons on a circle and do have kinetic term 
for the quarks, these could be Neumann boundaries as well. We hope to come back to these 
interesting questions in close future. 



6 Discussion 

In our analysis we have seen the following remarkable fact. The general framework of world-sheet 
RG analysis of the matrix quantum mechanics on a circle is seen to be capable of reproducing 
the critical exponent and the expected scaling laws of the c = 1 fixed point and the T-duality 
property respected by its singlet sector. Moreover, it indirectly reflects the presence of the non- 
singlet sector from the tendency of the change of the world-sheet free energy to change sign at a 
critical radius. It also indicates the presence of the boundaries in the 2D string theory and gives 
a flavor of capturing their dynamics. All these serve to be an initial understanding and progress 
towards our main goal, to develop a scheme for dealing with the non-singlet sectors (or the 
world-sheet vortices) which govern the physics of 2D string theory on small circle. In a separate 
publication we will describe how the behavior of the non-singlet sector can be explicitly 
captured in our RG scheme with a gauged matrix model, by introducing an appropriate gauge 
breaking term and tuning the fugacity of the vortices coupled to it (along with the usual tuning 
of the other parameters in the double scaling limit). In this case, the RG flows will turn out to 
be capable to explore the interesting fixed points beyond the c = 1 universality class. 

As we have observed, by tuning the parameters M 2 and g we arrive at RG trajectories with 
appropriate scaling of the bulk and the boundary cosmological constants, suggesting the presence 
of the various boundaries of the underlying 2D theory. To explicitly see these boundaries for 
each of the Dirichlet and Neumann cases one needs to study the behavior of the operator 
creating boundaries of finite length. The most relevant operator that controls this length is the 
boundary cosmological constant. Such finite length boundary operators are computed from the 
free fermion wave function. Using our RG scheme one can study the behavior of these operators 
for the c = 1 fixed point at any radius and try to see the all possible boundaries. The results 
at large R can be compared with the free fermion wave functions as a test. The study can be 
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further extended to include the effects of the non-singlets. We have work in progress in this 
area. 

We would now comment on the subtle issues on dealing the nonlocal terms in the action, 
in computing the determinant emerging from the adiabatic integration over the vector row and 
column of the original matrix variable. In our method, we are expanding the determinant in the 
Fourier space about <p(t) = (small field approximation). We tackle the non-local integrals by 
introducing 'center of mass' (large) and 'relative' (small) coordinates, expanding all the functions 
around the latter and then integrating over the 'relative' coordinates. We keep the contribution 
from the non- localities in the action up to the kinetic term n 2 (t), and the contribution from the 
higher order terms in (j)(t) up to the <j) 3 (t) term. The kinetic term is the most important nonlocal 
term to capture the effect of the world sheet vortices. As by the parametrization ( p.7|) it also 
gives rise to the kinetic term of the vector fields, it plays an important role in understanding 
the nonlocal Neumann boundaries, the Dl brane. The higher order terms like Tr(n m (t)0 n (t)) 
generated by expanding the determinant give rise to powers of traces and are redundant for 
the RG analysis. The higher order terms are also negligible due to small field approximation. 
However, it will be interesting to study the scheme with more couplings to see the contribution 
from the higher order terms and the convergence of the scheme. Also it is an interesting question 
to understand the relevant operators of the flow from the matrix quantum mechanics point of 
view. 

An interesting solution of the two-dimensional theory, apart from the flat space background 
with linear dilaton, is the two-dimensional black hole [^, ^TJ. It is partially believed that 
the nonperturbative formulation of two-dimensional string theory in terms of an integrable 
theory of noninteracting nonrelativistic fermions of the matrix quantum mechanics can deal 
issues like black hole evaporation and gravitational collapse. In any case, the condensation 
of world-sheet vortices of the two-dimensional noncritical string theory should describe two- 
dimensional black hole background, which had been a long standing challenge to obtain from 
matrix model. The recent work by Kazakov, Kostov and Kutasov gives proposals for a matrix 
model for two-dimensional black hole |6^| . In our forthcoming paper |63| , we will address this 
issue in more detail. Also in this context, one can consider integrating out several vector rows 
and columns simultaneously to see the effect of inserting multiple boundaries. This is analogous 
to the situation of the decay of several unstable branes and is anticipated by some authors to 



give black- hole like phase at the end point of the decay []45j , |64}| . However, the multiple row and 
column integration presumably would insert redundant operators of the form (Tr0^) n , which 
are not likely to give rise to black hole fixed point. 
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Appendix 

A The Feynman Diagrams 



Here we evaluate and discuss the terms in different orders of the series E[g, M, 4>n,R,N] fl2.21|) 
using the summation rules discussed in section (3.3) and the relation ( |2.13| ) for the inverse 
Fourier Transform. 



A.l The terms of order O(0) 



g Tr 



V 1 



+ M 2 ) 



2M 



2-kR 



coth(nMR) / dt Tr <j>(t) 



(A.l) 



A. 2 The terms of order O(00) 



9 -Tr 



E 



+ M2 ) (w + M2 ) 



g 2 f dt\dt2 



(2ttR¥ 



Tr (<P(h)<P(t 2 )) 



E 



m,n 



exp (i(n - m) h/R) exp (i(m - n) 



(A.2) 



Now, changing the variables to 'center of mass' and 'relative' coordinates defined respectively 
by 



T 



we have, 



h + t 2 



dti dt 2 = 2 ) dT dr = 2 dT dr . 



T 



h-t 2 



(A3) 
(A.4) 
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Hence, 



where, 



and, 



Tr 



E 



■'m—nYn—m 



2wR 



4tt 2 /? 2 J 



i-rrR 

dT / dr Tr(0(T + r)0(T - r)) ^ 



exp (i( 



n — m 



2t/R) 



Att 2 R 2 



2-kR 



ttR 
ttR 



dT / dr Tr(0(T) 2 - r 2 0(T) 2 ) ^ 

J-nR m „ 



(w + !)(& + 1) 
exp — to) 2r/i?) 
(f + l)(f + 1) 



g 2 F g2 (R,M) dTTr(<P(T) 2 /2) + g 2 F g2 (R,M) dT Tr(0(T) 2 /2) , (A.5) 
Jo Jo 



F ff2 (i?, M) 



7T_R 



exp (i(m — n) 2r/Rj 



2 - 2i?2 ^ ^(f + M2 )(f + M2 ) 



2tt_R 



2tT 2 j R 2 Jo 
1 



exp (i( 



m — n) t 



>/R) 



(r' = 2r) 



2M 2 sinh 2 vrMi? 



1 



M 3 sinh 2 ttM/? 



/ dr' cosh M(irR + r') cosh M(vri? - r') 
Jo 

-it MR cosh 2-nMR + - sinh 47rM.R 
2 8 



(A.6) 



F g2 (R,M) 



^ pttR 



exp (i(m -n) 2t/R) 



1 



2wR 



8n 2 R 2 

1 



exp (i( 



to — n r 



'/R) 



m 
r2-wR 



tiw + ^iw + M 2 )' 



(t> = 2r) 



8M 2 sinh 2 nMR 
1 

M 5 sinh 2 vrMi? 



/ dr' r' 2 cosh M(7ri? + r') cosh M^i? - r') 
Jo 

1 tc^M^R^ 
- — (1 + 8tt 2 M 2 R 2 ) sinhAnMR cosh27rM J R 

ttMR , , 

H cosh 47rMi? 

16 



(A.7) 
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A. 3 The terms of order O(000) 



9 -Tr 
6 



E 



^(* +M2 )(* +M2 )(* +M2 ) 



'm—n'rn—k'rk—m 



exp (i{n — m) t\/R) exp (i(k — n) t 2 /R) exp (i{m — k) t\j ' R) 



E 



m,n 



(A.8) 



Using redefinition of the variables into the "center of mass" and the "relative coordinates", 

T= ^(h+t 2 +t 3 ), T 1 = (t 1 -t 2 ), T 2 = {t 1 -t 3 ), 



dti dt 2 dt 3 = jf 2 ' dT dri dr 2 = dT dri dr 2 . 

\T,T U T 2 J 

Considering T\ and r 2 to be small and keeping the order O(0 3 ) term, above series could be 
evaluated as, 

,3 r 



6 







rJf + M2 )(i + M2 )(§+ M2 ) 

p2nR pttR pttR 

= iT LJ Tl LJ T2 



E 

m,n,k 



exp (inTi/ R) exp ( — imr 2 /R) exp (ik{r 2 — Ti)/ R) 



(f + M*)(f + M^)(f + M^) 

g 3 F g3 (R,M) / rfT Tr(0(T) 3 /3) 
./o 



(A.9) 
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where, 

F 9 s(g,R) 

exp (inTi/R) exp ( — imr 2 /Rj exp {ik{j2 — T\)/R) 



=wLj ti L dT2 



E 



1 



dn / dr 2 cosh(7rM J R - Mt x ) cosh(vrMi? + Mr 2 ) 

tt_R J-ttR 



16M 2 7ri?sinh 2 vrMi? J_ 
exp (ik(r 2 -n)/R) 



df I df ( cosh(27rMi? + Mr) + cosh Mf ) V . , ., 



(f + ^ 2 ) 



7lR ^ R 



df 



( cosh(27rMi? + Mr) + cosh Mr) cosh(7rMi? - Mf) 



16M 3 sinh 3 vrMi? Jo 
+ ( cosh(27r Mi? - Mf) + cosh Mf) cosh(7r Mi? + Mf ) 

7rMi? 

= = [47rMi? (3 cosh vrMi? + 2 cosh 3vrMi? 

64M 5 sinh 3 7rMi?(cosh 27rMi? + cosh 47rMi?) L v 

+2 cosh 57rMi? + cosh 7nMR) + sinh TrMi? + sinh 3vrMi? + sinh 5vrMi? + 2 sinh 7vrMi? 

+ sinh97rMi?] (A. 10) 

B The Scaling Dimensions 

Here we will summarize the general expession of the matrix Qk,i = ' an d its eigenvalues 
Qi and Q 2 , the scaling dimensions of the relevant operators, at different fixed points. 

l(5h-l) + 3g* 2 (F; 3 -^ 



d 
dM 

For the gaussian fixed point 



fi 11 = -(5/ i -i) + 3^(F; 3 --F; 2 ) 

12 y V dM 2 dM J 

n 21 = <?*((i/m* - m*)f; 2 - m*f; 2 ) 

n 22 = h + g* 2 ^ ((1/2M* - m*/2)f; 2 - M*/2i>) (B.ll) 



Q n = (5h- 1)/2 = -1/2, tt 22 = h = 0, tt 12 = tt 21 = 0. (B.12) 
For the nontrivial fixed point 

Q u = l + 5h/2 = 1, tt 22 ~h~0, f] 12 ~ft 21 ~0. (B.13) 
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